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Abstract 

In the present work the sensitivity of calculated /3/3-decay amplitudes to a realistic residual 



interaction is analyzed in the framework of the approach of Refs. [ij, ll3|]. Both the Gamow- 
Teller (GT) and Fermi (F) matrix elements M'^'^ for two- neutrino /3/3 decay (2i//3/3decay), along 
with the monopole transition contributions to the total matrix elements M^'^ of neutrinoless /3/3 
decay (Oi//3/3decay), are calculated within the quasiparticle random-phase approximation (QRPA). 



Decompositions of M^'^ and M^'^ are obtained by the method of Refs. Ij, [iSj] in terms of the 
corresponding energy-weighted sum rules S. It is shown that in most of the cases almost the 
whole dependence of M^*^ and M^'^ on the particle-particle (p-p) renormalization parameter gpp 
is accounted for by the gpp dependence of the corresponding sum rules S. General expressions 
relating S to a realistic residual particle-particle interaction are derived, which show a pronounced 
sensitivity of S to the singlet-channel interaction in the case of F transitions, and to the triplet- 
channel interaction in the case of GT transitions. Thus, the sensitivity of M^*^ and M^'^ to the 
SU(4)-symmetry-breaking part of the p-p residual interaction is dictated by the generic structure 
of the /3/3-decay amplitudes. Therefore, a choice of this part in a particular calculation needs a 
special caution. Finally, a better isospin-consistent way of renormalization of a realistic residual 
p-p interaction to use in QRPA calculations is suggested. 
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I. INTRODUCTION 



Study of neutrinoless double beta (Oi//3/3) decay provides one of a few ways to probe the 
absolute neutrino mass scale with a high sensitivity [H, H]- To be able to deduce the effective 
Majorana neutrino mass from the measured half-lives of the decay, reliably calculated nuclear 
matrix elements (NME) M"*^ are needed in addition. 

The NME have been calculated in different approaches: the quasiparticle random- 
phase approximation (QRPA) [3h5|] (including the recent version of the QRPA accounting for 
deformation [^), the nuclear shell model (SM) 0], the projected Hartree-Fock-Bogoliubov 
method (PHFB) ^ , the interacting boson model (IBM-2) [q!] , and the generator coordinate 
method with particle number and angular momentum projection (GCM-I-PNAMP) [lo(]. 
There has been great progress in the calculations over the last decade, and now M"*^ of 
different groups (apart from the SM one) seem to converge. However, M'^" of the SM are 
systematically and substantially (up to a factor of 2) smaller than the results of the other 
approaches. This discrepancy calls for a better understanding of the aspects of nuclear 
structure that affect largely the calculated NME. 

All the QRPA calculations of the nuclear /3/3-decay amplitudes have revealed their sensi- 
tivity to the strength of the particle-particle (p-p) interaction in the triplet channel,^ which 
is mainly due to a high Qpp sensitivity of the contribution of transitions through the inter- 
mediate l"*" states. Such a behavior of the two- neutrino 2i//3/3-decay NME M'^'^ was found 
for the first time in Ref. Uj. Reference [llj contains a clue that the sensitivity may be 



related to the restoration of the Wigner spin-isospin SU(4) symmetry in nuclei. The authors 
of Ref. [HI also showed that the sensitivity is not an artifact of the QRPA but also shows 
up in an exactly soluble schematic model. 

An idea to use the concept of softly broken SU(4) symmetry as a basis for describing 
2u(3(3-decaj amplitude was put forward in Ref. [l^. Because M^*^ vanishes in the limit 
when the SU(4) symmetry is exact, it is natural to express M"^" explicitly in terms of those 
parts of the nuclear Hamiltonian H that are responsible for the violation of the symmetry. 
An identity transformation introduced in Ref. 13] allows one to shed light on the general 



properties of the 2i//3/3-decay amplitude. However, the computational realization of this idea 
in Ref. [12] made use of a oversimplified model of independent quasiparticles and, therefore, 
could not address the question of th e q„ p sensitivity. 

The next step was made in Ref. where the basic concept of Ref. [l2| was applied 



in the framework of a QRPA model. As in Ref. the starting point for the analysis 
was a model-independent, identity, transformation of M^'^ . That allowed the authors to 
partition M^'" into two terms which are sensitive to different parts of H. The dominating 



source of the gpp sensitivity was associated in Ref. [13| with a specific energy-weighted sum 
rule 5" for double beta decay, that depends exclusivel y o n the residual p-p interaction. An 
analytical representation for S was obtained in Ref. 1^ within the QRPA for the case of 
a simple separable p-p interaction, which showed that S vanishes at the point where the 
SU(4) symmetry is restored in the p-p sector of the model Hamiltonian. The rest of M^^ 
was shown to behave smoother on Qpp for the realistic values of the p-p interaction strength, 
and was mainly determined by the other important source of breaking the SU(4) symmetry, 
namely the spin-orbit part of the nuclear mean field. 



^ The so-called "gpp problem" , whereby the factor gpp is used to renormalize the p-p interaction in the 
QRPA equations 
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In the present work we apply the approach of Refs. [12|, Il3| to analyze the sensitivity of 
calculated /3/3-decay amplitudes to a realistic residual interaction. Both the Gamow- Teller 
(GT) and Fermi (F) matrix elements M'^^ , along with the monopole transition contributions 
to M°^, are calculated within the QRPA, making use of a realistic residual interaction (the 
Brueckner G matrix) as described in Refs. 0-[5[. General expressions relating the sum rule 
S" to a realistic residual p-p interaction are derived, which show a pronounced sensitivity of 
S to the singlet-channel of a two-body interaction in the case of F transitions, and to the 
triplet-channel of a two-body interaction in the case of GT transitions. In this connection, 
S for GT transitions would be a better quantity for fitting Qpp than the experimental M^^, 
provided S could be measured. ^. Identity partitions of M^'^, as well as the monopole 
transition contributions to M^'^, are obtained by the method of Refs. (12I . 13| . It is shown 
that in most of the cases almost the whole gpp dependence of M'^^ and M^'^ can be attributed 
to the Qpp dependence of the corresponding sum rules S. Thus, the sensitivity of M"^^ and 
M^'^ to the SU(4)-symmetry-breaking part of the p-p residual interaction is unavoidable 
since it is simply dictated by the generic structure of the (3(3 amplitudes. Finally, a better 
isospin-consistent way of a renormalization of a realistic residual p-p interaction to use in 
QRPA calculations is suggested. 



II. IDENTITY PARTITION OF /3/3-DECAY AMPLITUDES 



We adopt here the same line of reasoning as presented in Refs. |l2l . |l3|. The F and 
GT transitions are treated in a uniform way, and the single-particle (s.p.) operator 
/3j = governs allowed Fermi {g,j=Q = 1) or GT {gj=i = cr) /3-transitions, 

respectively. Here, J = 0, 1 is the angular momentum of a state (with positive parity) that 
can be connected by the operator /3j with the ground state of an even-even nucleus. 

The 2z//3~/3~-decay amplitude can be written in the form [2| 

Mf = Y^^ . (1) 

s 

Here, gs = {0'f\\(3j\\J~^, s) and fs = {J~^ , s\\(3j\\0f) are the one-leg transition matrix elements 

of the operator P~ between an intermediate state s of the isobaric nucleus (A^ — 1,Z + 1) 
and the ground states (g.s.) of the parent {Z, N) and the final {Z + 2, N — 2) nucleus, 
respectively. The energy denominator Cjs in Eq. ([T]) is the excitation energy of the sth state 
relative to the mean g.s. energy of the initial and final nucleus, Us = {cOs{i) + i^s{/))/2, with 
Us(i) = Eg — Eq. {uJs{f) = Eg — Eq^) representing the excitation energy of the s'th state 
relative to the g.s. of the initial (final) nucleus. 

The following partition of M^*^ ([1]) can be performed [IJ 



Mf = + % (2) 



ri2 
^9 



^ Realistically, the absolute value of S can only be determined experimentally (from charge-exchange reac- 
tions or single-/? decays) if the single-state dominance is realized in 2j//3/3-decay of one or another nuclear 
system. 
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Sf 



E 



(3) 
(4) 



where Qg is, in principle, an arbitrary energy. Note, that the second term in the r.h.s. of 
Eq. ([2]) is proportional to S'^ , which has the form (j4]) of an energy- weighted sum rule. One 
can get the following expression for S"] 



'Til 



^r = -^(0;|5j-|0+) , SJ 



/3,7, A7,^ 



(5) 



where H is the nuclear Hamiltonian and the tilde denotes the time-reversal operation. A sum 
over all the spin components is assumed in the case of the GT transitions. The expression 
([5]) contains explicit information about symmetry properties of the nuclear Hamiltonian in 
terms of the corresponding commutators. If one chooses ujg in Eq. ([2]) to coincide with 
the energy log of the corresponding giant resonance (the isobaric analog state (IAS) or the 
Gamow- Teller resonance (GTR)), one sees immediately that both M''^j and S'f' vanish in 
the limit of the exact isospin SU(2) (J = 0) or spin-isospin SU(4) (J = 1) symmetries. It can 
be seen for small deviations from a symmetry that S']'^ depends linearly on the symmetry- 
breaking terms of the Hamiltonian, whereas M''^J' has a weaker, quadratic, dependence on 
them. 



It is shown in Ref. [13| that M'j and Sf are sensitive to different terms in H. By making 
use of the quasiboson approximation (QBA), S"]" was demonstrated to be only determined 
by the p-p part of the residual interaction (chosen in that work in a separable form). An 
analytical expression for the GT sum rule Sl^ was also derived in Ref. 13| : 



St 



2u 



3^^(l-(? 



Go 



(6) 



where A„ and Ap are the pairing gaps for neutrons and protons, and g'^ = ^ is a ratio 
of the strength of the triplet (spin 5 = 1) p-p interaction G\ and the singlet (spin 5 = 0) 
one Go (the latter governs the pairing correlations in nuclei). In the derivation of Eq. ([6]) 
the BCS vacua were taken the same for initial and final nuclei. The point g'^^ = 1, where 
Sf'^ = 0, corresponds to the restoration of the SU(4) symmetry in the p-p sector of the 



model Hamiltonian of Ref. 13 



It is useful to consider also a closure 2z//3/3-decay matrix element M'^'^ 



J {cl} 



^^■^j {cl} 



(7) 



and to apply to it a partition similar to Eq. 



{cl} 



M"f{.i} + 



M' 



J {cl} 



E 



{ujg - UJs)gsfs 



(8) 
(9) 
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The nuclear matrix element of Oz//3/3decay is given by a sum of the partial amplitudes 
M^'^{J) of transitions via the intermediate states of all multipolarities J'" (see, e.g., Refs. 0- 

9): 

M"^" = M^^iJ)- M'^'iJ) = J2 M^'iJ) (10) 

J s 

Mj" also can be partitioned: 

MO-(J) = M'°"(J) + ^ (11) 

Ug 

^^ {Ug - UJ,)M^/{J) _ SOU ^J2-^^0u^j^ (12) 

s s 

Below only the contributions M^'^{0~^), M°'^(l^) of the monopole transitions through the 
intermediate states O"*", l"*" are analyzed, which are known to be sensitive to the p-p interac- 
tion. 

We stress again that the all the transformations of the /3/3-decay amplitudes Eqs.(p|)- (fT2|) 
introduced above are identical, and therefore do not rely on any nuclear model. 



III. DERIVATION OF Sf FOR A REALISTIC RESIDUAL INTERACTION 

In this section we present a derivation of 5*7" and ([5]) in the case of a general realistic 
residual interaction. 

One sees immediately that a single-particle mean field, as containing only isoscalar and 
isovector terms, exactly drops out of the double commutator defining Sj~ . Only the 
residual two-body interaction V = | ^ fab contributes (we work here in the first quantiza- 

tion since it simplifies the further derivation): 
SJ- = iT^Sfiab) (13) 

Sf{ab) = [f3-f,^j, Vab\] = f3ab,jf^ab,jVab + VabP'b/Pab^J " Pab,jVabPab,J - Pab,jVabPab,J i'^^) 

Here, /3± ^ = gj{a)T^{a)+gj{h)T^{h) = gjT^ + gJt^, gj = ^{gjia)+gj{b)), gj = ^{gAa)- 
gj{b)), T± = r±(a) + r±(6), t± = r±(a) - r±(6), and = r,(a) + r,(6). 

Thus, the problem of calculating Sj~ ([5]) is reduced to the problem of calculating the 
two-body double commutator Sj'^{ab) f|T^ . Bearing in mind a further calculation of Sj'^ as 
a matrix element of Sj~ between antisymmetric nuclear wave functions of the initial and 
final nuclei, the operators of different permutation symmetry gf and T^, acting in the 
two-body space are introduced. 

By making use of isospin projection operators, the original two-body interaction, which 
is considered exactly isospin symmetric here, can be partitioned into two components, cor- 
responding to different projections of the total isospin of the two-nucleon system: 

Vab = Vab{Tz = 0) + Vab{\Tz\ = 1) (15) 
Vab{Tz = 0) = Vab , Vab{\Tz\ = 1) = Vab ' , (16) 
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each of which can further be represented in terms of the interaction components correspond- 
ing to a definite total isospin T = 0, 1: VahiT^ = 0) = VabiT = l,Tj = 0)nT=i + Vab{T = 
0,Tz = 0)11^=0; Vab{\Tz\ = 1) = Vab(T = 1, \Tz\ = 1) (Ilr is a corresponding total isospin 
projection operator). 

The standard way of renormalization of the residual p-p interaction in the proton-neutron 
QRPA is: 

Vab gppVab{Tz = 0)+ gpairVabi\Tz\ = 1), (17) 

because only Vab{\Tz\ = 1) enters the BCS gap equations, while VabiT^ = 0) is responsible 
for the mixing of proton-neutron excitations (here for the sake of simplicity we consider the 
same Qpair for proton and neutron subsystems). 

For the F transitions {go{a) = 1) one arrives at the following expression (see Appendix): 

So{ab) = 2{gpair - Qpp) (T") Vab{T = 1, = 0). (18) 

This result shows that the renormalization f lT7|) obviously breaks the original isospin symme- 
try of the residual interaction if gpp ^ Qpair (that is usually the case in most of the realistic 
QRPA calculations; we shall later how this drawback can easily be remedied by a different 
renormalization that is more isospin consistent than the one of Eq. ( !T7|) ). 

The corresponding expression for the GT transitions {gi = cr) is more involved (see 
Appendix) and reads 

Sf'iab, S = 0) + Sl^iab, S = 1) (19) 

6 (r-)' [gparrVabiT = 1, ^ = 0) - gppVab{T = 0, S = 1)] Us=o (20) 
2 (T^^ [{gpp - gpair)Vab{T = 1, = 1) 

+gpp{vab{T = 1, 5 = 1) - Vab{T = 0, 5 = 0))] (1 - ns=o), (21) 

where the operator 115=0 = |00)(00|g projects onto the spin S* = state of two nucleons. 
The operators Vab(T, S) = {TS\vab\TS) are the expectation values of the original two-body 
interaction in two-body spin-isospin states \TS), and therefore only depend on the spatial 
coordinates of two nucleons. 

Now we proceed with calculations of Sj'^ as the matrix elements of the two-body operators 
Sj^ between the ground states of the initial and final nuclei. Since Vab is of a short range, 
then the nucleon pairs in the relative spatial s wave must predominantly contribute to the 
Sj^. The operator (T~)'^ in Eqs. (ITSl and [T9|) transforms a T = 1 neutron pair into a T = 1 
proton pair [and becomes just a number (T~)^ = 2 for these states]. Such pairs of nucleons 
must then be in the state with the total spin 5* = to assure antisymmetry of the total two- 
body wave function. This means that the term Si" {ah, S = 1) fl^ as projecting onto S = 1 
states can safely be neglected. Further, one can anticipate in advance that the dominating 
contribution to Sf" should come from the paired neutrons and protons in the two-body 
state J'^ = 0"^. By taking into account only this leading contribution of paired nucleons, one 
arrives at the following representation for S'^'^: 



Sl''{ab,S = 0) = 
Sf''{ab,S= 1) = 



^2u /I f \ nipair) 

'0 — U ~ 9pp)9pairOQ 

^(pa^r) ^ ^ J = Q, T = 1) (0+ 1 [cjcj] OO [c„C JoO 1 0+) , (23) 



'S'o'' — (1 9pp)9pairSQ^ , (22) 
(pair) 



pn 
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where g'^^ = gpp/gpair, and [cJ4]oo and [c„c„]oo are bifermionic operators made of the coupled 
to J'" = 0^ particle creation and annihilation operators c|,Ct {t = p,n), and the G matrix 
G{J = 0,T = 1) corresponds to the two-body interaction Vab{T = 1) (from which only the 
component Vab{T = 1, = 0) is active in the 0"*" channel). 
The corresponding expression for Sf'^ reads 

Sr = {l-lM9pa^rS?-\ (24) 
gipair) _ ^gipair) ^ _ gipp) ^ gipair) ^25') 

Si""'^ = EpnG'ippnn; J = 0,T = 1)(0+| [ctct]oo[cnC„]oo|0+), (26) 

where the G matrix G'{J = 0,T = 1) corresponds to the two-body interaction Vab{T = 
0, 5 = l)n5=onT=i in Eq.([20]). 3 

Let us note that retaining only the contribution of J = paired nucleons in the g.s. 
wave functions corresponds to the calculation of the double commutator ([5]), defining 5'^'^, 
in the QBA. In the QRPA this means that the commutator of two bifermionic operators is 
substituted by its expectation value in the BCS state, which is a c-number (cf. a derivation 
of the QRPA matrices A and B, e.g., Ref. [l8|). In fact, there are additional, beyond the 
QRPA, contributions to 5"^^ from the pairs with J > in the correlated g.s., but they must 
be suppressed as the following arguments suggest. For the pairs with J > 0, which are 
in the relative s wave and may therefore essentially contribute to the sum rule, the total 
J must then coincide with the total orbital momentum of the pair. However, as previous 
calculations of the Oz//3/3-decay transition densities have shown 0, the contributions 



from J > pairs come from larger internucleon distances, substantially exceeding the short 
range of the A^A^-potential, then that of J = (the latter peaks essentially at 1-2 fm). This 
in combination with the short range of Vab causes a suppression of the J > contributions. 

An estimate of such a suppression is beyond the scope of this paper and deserves a 
separate study. These contributions can be estimated if in a calculation of Oi//3/3-decay one 
substitutes the Coulomb-like r-dependence of the neutrino potential by a corresponding r- 
dependence of the T = 1 component of a A^A^ potential. Note, that even in the case of the 
long-range neutrino potential the 0"'"-pair contribution is by far the largest one. 

For the BCS description of pairing, by taking the same BCS solution for the initial 
and final nuclei, one gets (Oj | [cJ4]oo|0^) = (0^ | [ciCi]oo|0^) = jtUtVt, where u,v are the 
Bogoliubov coefficients, UtVt = and the pairing gaps satisfy the gap equation: 

A, = ^ ^ Gittt't'; J = 0, T = l)hu,v,. (27) 
Jt t' 

It is noteworthy that the same expressions fl22p - (l2^ for can be obtained by considering 
the double commutator ([5]) in the QBA, like it was used in the derivation of Ref. [l3|. A 
compacter expression can be obtained for the pairing sum rule s'f"'^'^^ by taking into account 
the gap equation Eq. ( 1271) : 

9pa.r Sr-^ = \ T.Wj\\n? ^^^"(^^+-^") (28) 

^ pn ^p^n 



^ In such a deuteron-type channel with T = 0, S* = 1, because of the tensor mteraction, the attraction 
must be stronger than in the one with T — 1,5 = 0, thus one can expect 7^ > 1 in Eq. (j25p prior to any 
calculation. 
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IV. RESULTS AND ANALYSIS 



Within the QRPA approach of Refs. jl-lsl we have computed Mj" (P), 5f ©, Mj^^^^y © 
along with M^^i-J) and ^°^(J) for the decays ^^Ge^^Se, iooMo^i°°Ru, and 
^^°Te— J-^^^Xe ( J = and J = 1 for the F and GT transitions, respectively). The parametriza- 
tion of the Woods-Saxon mean field is adopted from the spherical calculations of Refs. 
For each of the nuclei two sizes of the s.p. basis, the small one (s.b.) and the large one (l.b.) 
in the notation of Refs. [sl-lsj, are used in the calculations. The small s.p. space consists of 9 
levels (oscillator shells A^=3,4) for A = 76 and 13 levels (oscillator shells A^=3,4 plus f + h 
states from = 5) for A = 100, 130. The large basis contains 21 levels for A = 76, 100 
(all states from shells = 1 — 5), and 23 levels for A = 130 (A^ = 1 — 5 and i orbits from 
N = 6). 

As in Refs. [3-:5'3; the nuclear Brueckner G matrix, a solution of the Bethe-Goldstone 
equation with different nucleon-nucleon potentials (Bonn-CD, Bonn-C, Argonne V18 and 
Nijmegen I), is used as a residual two-body interaction. The results obtained with different 
G matrices look pretty similar to each other, and in all the figures below only the results 
obtained with the Bonn-CD G matrix are represented. 

First, the BCS equations are solved to obtain the Bogoliubov coefficients u and v, the 
pairing gaps A, and the chemical potentials. To correctly reproduce the experimental odd- 
even nuclear mass differences for both protons and neutrons in initial and final nuclei, four 
slightly different renormalization factors Qpair in Eq. (!27|) are needed. Here, we approximate 
the single parameter gpair of the preceding section by taking the average value of these four 
factors. 

As in Refs. we set the particle-hole renormalization factor Qph = 1 in the QRPA 

equations. The calculated energy of the giant GT resonance, which is essentially independent 
of the size of the s.p. basis, is well reproduced with such a choice of Qph [S-ll]. One must say 
that a particular choice of gph in the QRPA has no effect on the sum rules Sf'^ and Sq'^ (jlj), 
and S^^IGT) and S^'^{F) f|T2|) . This can be seen from the general analytic expressions fl22|24p 
(determined exclusively by the p-p interaction) and is confirmed by the direct calculations. 
The sum rules Sf (gD, and 5°^ (J) ([12]) calculated within the QRPA approach of Ref. 0- 
are shown in Fig. [TJ ^ One can see a universal character of the almost perfectly linear 
dependencies S^g'pp) for all the nuclei in question, with only a little dependence on the s.p. 
basis size. Only for ^^°Te— )-^^°Xe does one observe a slight depletion of S" as a result of 
weaker proton pairing in ^^''Te (only 2 protons upon the Z = 50 core). The fact that the 
dependence S{g'p^ is almost the same for different nuclei can easily be understood from 
the phenomenological A dependencies A oc and Go put into S oc AnAp/Go, 

Eq. ([6]), for separable forces. These estimates also explain why S{0) slightly increases with 
an enlargement of the s.p. basis [a smaller Qpair (or Go) is needed to fit the experimental 
A„, Ap]. 

Also, Sf' and the corresponding S^'^{J) have a surprisingly good quantitative agreement 
(a possible clarification of this numerical observation needs a further study that is beyond 
the scope of this work). All the Sp lines cross zero very close to g'^^ = 1, in an excellent 
agreement with Eq. fl22|) , despite the aforementioned numerical differences in the fitted gpair 



* All the Oi^ quantities are calculated here by retaining only the leading Coulomb-like radial dependence in 
the neutrino potential 
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FIG. 1: (Color online) The GT and F sum-rules Sf" and S^" (gD (the solid and dashed lines in 
left column, respectively), and £''''"(1) and S^'^{0) (fT^ (the solid and dashed lines in right column, 
respectively) calculated within the QRPA approach of Ref. 0-0] • The results for the small s.p. 
basis (s.b.) are represented by the red lines, and for the large one (l.b.) are represented by the 
black lines. 



for protons and neutrons and the overlap factor used in the present QRPA calculation. The 
ratio Si = SSp at g'^^ = comes out almost exact in the calculation, again in excellent 
accord with Eq. (!25l) . All the Si lines cross zero at g'pp = < 1, as an expected result of 
a stronger attraction in the deutron-type channel T = 0, S = 1. The calculated values of 7^ 
are listed in Table [H here for different choices of the residual interaction. ^ For a given choice 
of the residual interaction, one sees again an impressive universality of this parameter, with 
the largest value 7^ ?a 1.5 in the case of the Bonn-CD A^A^-potential and the smallest value 
7j 1.3 in the case of the Bonn-C A^A^-potential. 

In this connection, provided 5* for GT transitions could be measured, it, as depending 
exclusively on the residual p-p interaction, would be a better quantity for fitting gpp than the 
experimental M2j/-ea;p fact, the absolute value of S can be determined experimentally 
from charge-exchange reactions or single-/? decays if the single-state dominance for 2i//3/3- 
decay is realized in one or another intermediate nuclear system. 

In Fig. [2] we plot Mj'^^g'pp) (solid line), S^/u"^ (dashed line), and their difference, 
M''^j{g'pp) (dot-dashed line), calculated within the QRPA according to Eqs. The 
upper and lower panels contain the results for F and GT transitions, respectively. The value 
of ujg is calculated as the mean energy of the GT or F strength distribution in the /3~ channel 
with a low-energy cut-off of 10 MeV. One observes that the dependencies Mj^i^g' ) for differ- 



^ Here the same BCS vacuum of the initial nucleus is used also for the final one to avoid the influence of 
the overlap factor. The latter spoils the observed universality of 7^ in Table U by a few percent 
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TABLE I: The calculated parameter 7^ Eq. (|25p for different nuclei, basis sizes and choices of the 
residual interaction. 





Bonn CD 


Argonne VI 8 


Nijmegen I 


Bonn C 


™Ge (l.b.) 


1.485 


1.418 


1.371 


1.274 


™Ge (s.b.) 


1.477 


1.413 


1.368 


1.280 


iooMo (l.b.) 


1.493 


1.426 


1.378 


1.279 


ioomo (s.b.) 


1.498 


1.426 


1.378 


1.280 


i30Te (l.b.) 


1.485 


1.419 


1.371 


1.276 


i30Te (s.b.) 


1.481 


1.416 


1.370 


1.280 



ent basis sizes, s.b. and Lb., plotted as functions of g'^^ look much more similar to each other 
than in the usual case when they are represented as functions of gpp. As for the function 
S/Ug, it shows basically a linear dependence on g'pp governed by the corresponding behavior 
of S. One can see from the figure that the dependence M''^f{gpp) is much smoother than the 
original one, Mj'^^g'pp), for realistic values of g'pp (apart from one exception — Mf'^^g'pp) for 
^°°Mo-)-^°°Ru, where one approaches very close to the point of the QRPA collapse). Also, 
all Mq'^ cross zero very close to the point g'pp = 1 of restoration of the isospin symmetry of 
the renormalized residual interaction f lT7|) . 

In Fig. |3]we show for completeness the closure matrix elements Mj^^^i^, Mf" (solid lines), 

and Sf' /ug, Sj^/ug (dashed lines) and the differences, M'"^/ and M'^f (dot-dashed lines), 
calculated within the QRPA according to Eqs. dB-® (2z//3/3) and Eqs. ([II])-(II2D (Oz//3/3) in 
s.b.. The upper and lower panels again contain the results for F and GT transitions, respec- 
tively, and in each panel the left and right columns of figures show the results for the 2z//3/9 
and Oz//3/3 cases, respectively. One sees that in this case almost the whole g'pp dependence of 
^T{ci} is governed by the g'pp dependence of and S'j', respectively. 

The thin solid horizontal lines in the lower panels of Fig. H] represent the corresponding 
experimental values M^''"^^^ obtained in Ref. [15|] with the unquenched value of the axial- 
vector coupling constant gA = 1.25. Here one immediately sees a problem — one needs 
g'pp < 1 to fit Mi"^'^^^, but for such a choice the renormalized residual interaction f lT7|) 
breaks isospin and M^'^ are spuriously large. However, this drawback is very easy to remedy 
by a slightly different prescription of a renormalization of the residual interaction. According 
to Eqs. fl22| 121]), the F transitions are sensitive to gppVab(T = 1, S = 0), whereas the GT 
transitions are sensitive to gppVab{T = 0,S = 1). Therefore, it suffices to renormalize different 
T components of Vab as 

Vab gppVab{T = 0) + gpairVab{T = 1), (29) 

to have an isospin symmetric interaction, which allows one at the same time to fit the 
odd-even nuclear mass differences by means of gpair and M^*^"*^^^ by means of gpp. Direct 
QRPA calculations using the renormalization Eq. f l2I?]) show that the GT g'pp dependencies 
shown in Figs. [TH3] stay practically the same (the change would be barely visible in the 
figures), whereas the the F g' dependencies become constant, equal to the corresponding 
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FIG. 2: (Color online) Mf{g'pp) (solid line); Sf/ul (dashed line); and their difference, M'j{g'pp) 
(dot-dashed line), calculated within the QRPA according to Eqs. ((l])-((4l) with different basis sizes, 
s.b. and l.b.. The upper and lower panels contain the results for F and GT transitions, respectively. 
The thin solid horizontal lines in the lower panels represent the experimental values M^''"^^^ 
obtained in Ref. with the unquenched value of the axial- vector coupling constant qa = 1.25. 
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values at g'^^ = 1 in Figs. [IH31 The effect of this new way of renormahzation of the residual 
interaction on the total Oi//3/3-decay NME will be investigated elsewhere, but one can 
already anticipate that Mp^ will come out slightly smaller than in Ref. jsj-jsl, whereas Mq'j, 
will barely be affected. 

The results discussed above demonstrate very little dependence on the s.p. basis size. 
However, the smallest QRPA s.p. basis used in the analysis, 2huj for A = 76, is still 
much larger than the corresponding one of the SM which contains only four s.p. levels: 
1^3/2) 0/5/2, lpi/2 and Ogg/2 [3]. The problem with this small Ohu SM basis is that the spin- 
orbit partners O/7/2 and OgT/2 are missing, which leads to a strong violation of the Ikeda sum 



rule (ISR) [16| (the QRPA satisfies the ISR exactly, see a detailed discussion in Ref. |17l|). 

It is instructive to see what happens to the sum rules S if one uses the SM s.p. basis. 
According to the arguments after Eq. fl2T|) . S are mostl y d etermined by the paired nucleons 



with J = in the g.s.. One can see in Fig. 5 of Ref. [17| that the contribution of J = 
pairs to M'^'^ when calculated in the QRPA with the SM s.p. basis comes out in fairly good 
agreement with the corresponding result of a genuine SM calculation. This fact makes us 
confident that a QRPA calculation of the sum rules S with the SM s.p. basis should give a 
reasonably good estimate for the corresponding SM result. 

The results of such a calculation are listed in Table [Tll It can be seen that the F sum rule 
So is still reproduced well in the SM basis. However, the GT sum rules Si"""''^ and ^^^^ come 
out strongly underestimated as a result of the missing contributions to 5*1 from the spin-flip 
transitions involving O/7/2 and OgT/2 s.p. states. The same happens to the parameter 7^, 
which is almost three times too small. 

Thus, one may state that the original inherent sensitivity of S to the SU(4)-breaking part 
of the residual p-p interaction gets spuriously weak in the SM basis for A = 76 system. To 
restore it, one must include the missing spin-orbit partners to the SM s.p. space, as, for 
instance, was the case in the SM description of decay of ^^Ca. The usual argument of 
the SM, that the weights of the admixtures of the missing states in the g.s. wave function 
are small, does not work here, because the relatively small weights get compensated by large 
transition matrix elements to those states while calculating 5*. 

TABLE II: Comparison of different S calculated within the QRPA for ^^Ge in the SM basis ("4 
levels") and s.b. (="9 levels") (the Argonne-V18 G matrix is used). 

4 levels 9 levels 

4^"'''^ 13.8 15.8 

S'f"''^^ 21.9 46.9 

S^f"^ 12.3 66.3 

7, 0.56 1.41 



To conclude the analysis of this section, a calculation of M^'^ within the renormalized 
QRPA (RQRPA) 0-[5[ is performed and compared with the corresponding QRPA calculation 
(Fig. H]). One can see that the calculated RQRPA dependence MQ^^g'^p) does not cross zero 
at the physical point g'^^ = 1 of the restoration of the isospin symmetry as the QRPA results 
do. Thus, this adds another drawback of the RQRPA to the well-known violation of the 
ISR. 



12 




> 



A=76, s.b. 


A=76, S.b. 


-1, 

■- m' > 

1,1 


1,1, 


A=100, s.b. 


A=100,s.b. 


1,1 


1,1, 


A=130, s.b. 


A=130,s.b. 


1,1 


1,1, 



2 


-2 
2 

' 

-2 

2 



-2 



0.5 



0.5 



FIG. 3: (Color online) The matrix elements Mf^^^^, M^^iJ) (solid lines), Sf/ujg, S^''{J)/ojg 

(dashed lines) and the difference, M'j^ M'^" {J) (dot-dashed lines), calculated within the 
QRPA according to Eqs. dZ])-® (2z/) and Eqs. ([II])-([l2]) {2v) in the small basis (s.b.). The upper 
and lower panels contain the results for Fermi and Gamow- Teller transitions, respectively, and in 
each panel left and right columns of figures show the results for the 2v and cases, respectively. 



13 




FIG. 4: (Color online) M^" calculated within the RQRPA in comparison with the corresponding 
QRPA calculation. 



V. CONCLUSION 

In the present work the sensitivity of the calculated /3/3-decay amplitudes to a realistic 



residual interaction is analyzed in the framework of the approach of Refs. [12|, Il3|. Both 
the GT and F matrix elements M"^" for 2z//3/3decay, along with the monopole transition 
contributions to the total matrix elements M^'^ of Oi//3/3decay, are calculated within the 

in 



QRPA. Decompositions of M'^'^ and are obtained by the method of Refs. [12, ll3 
terms of the corresponding energy-weighted sum rules S. It is shown that in most of the 
cases almost the whole dependence of M"^" and M^^ on the renormalization parameter Qpp is 
accounted for by the gpp dependence of the corresponding sum rules S. General expressions 
relating S* to a realistic residual p-p interaction are derived, which show a pronounced 
sensitivity of S to the singlet-channel interaction in the case of F transitions and to the 
triplet-channel interaction in the case of GT transitions. In this connection, S would provide 
the best quantity for fitting gpp if it could be measured (realistically, this can be done only 
in the case when the single-state dominance is realized in 2i//3/3decay). Thus, the sensitivity 
of M"^" and M^" to the SU(4)-symmetry-breaking part of the p-p residual interaction is 
dictated by the generic structure of the /3/3-decay amplitudes. Therefore, the choice of this 
part in a particular model and a further accurate calculation of its contribution to S needs 
special attention. Finally, a better isospin-consistent way of renormalization of the realistic 
residual particle-particle interaction to use in QRPA calculations is suggested. 
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Appendix A: Calculation of Sj{ab) 

In the isospin space of the two-nucleon system \TTz) {T = 0, 1; = 0, ±T), the isospin 
operators T~ and t~ can be expressed in terms of the isospin projection operators: 

T- = V2(|10)(ll|^ + |1 - 1)(10|^); r = V2(|00)(ll|^- |1- l)(00g (Al) 

It is easy to verify that T^Vabt^ = t^VabT^ = as a consequence of the isospin conserva- 
tion by strong interaction. In addition, one has Tz{a)Tz{h) = -^^^^^ — 1, Tz{a)Tz{h)T^ = 
— and t'Tz{a)Tz{h)t'' = (T^)^. Then one finds that T~~Vab{Tz = 0)T^ = Vab{T = 

1,T, = 0) (T-)^ t-Vab{Tz = 0)t- = -v,,{T = 0,T, = 0) (T-)^ and T-Vab{Tz = l)T~ = 
t~Vab{Tz = = 0, and finally arrives at Eq. f|T8l) for the F transitions. 

For the GT transitions gi{a) = cr^, and the operator gi{a)gi{b) = —cTaCrb has the eigenval- 
ues — 2(S'(S' + 1) — 3) in the spin space of the two-nucleon system \SM), S* = 0, 1; M = 0, ±5". 
The operators Qi and Qi are equal to S = \{,(^a + cTfc) and s = \{(Ta — (Tb), respectively , 
whose components can also be expressed in terms of spin projection operators: 

5~ = ^(|10)(ll|, + |l-l)(10y (A2) 

5+ = -^(|ll)(10|, + |10)(l-iy (A3) 
Sz = i2M\lM){lMl (A4) 

s 

and 

5" = ^(|00)(ll|,-|l-l)(00y (A5) 

s+ = i=(|ll)(00|,-|00)(l-iy (A6) 

s, = |10)(00|, + |00)(10|,. (A7) 
Then one finds the following expressions: 

g^Vab{T = 0)g^ = SzVabiT = 0)s, + 2{s-Vab{T = 0)s+ + s+Vab{T = 0)s-) (A8) 



Y,{lM\vab{T = 0)\lM)s 



\00){00\, + {00\vab{T = 0)\00)s 



J2\1M){1M\, 



(A9) 



and 



~gtvab{T = 0)gt = SzVab{T = 1)5, + 2{S~VabiT = 1)S+ + S^Vab{T = 1)S-) (AlO) 
= 2j2{lM\vab{T = 1)|1M)5|1M)(1M|, = 2vab{T =1,S = 1)(1 - |00)(00y. (All) 

s 
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Finally, having collected all the contributions, one arrives at Eq. fll9p for the GT transitions. 
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